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Abstract.  Transient  to  steady  motions  of  a  vapor  caused  by  the  evaporation  and  condensation  processes  occurring  at  the 
condensed  phases  placed  in  parallel  have  been  studied  based  on  the  Boltzmann  equation  of  BGK  type.  As  the  internal 
structures  of  the  condensed  phases,  the  temperature  fields  are  taken  into  account.  Because  of  this,  the  temperatures  of  the 
interfaces  become  unknown  parameters  and,  therefore,  the  condition  of  the  continuity  of  energy  flow  across  the  interface  has 
to  be  imposed  simultaneously  with  the  conditions  so  far  used  for  the  cases  with  no  internal  structures.  This  extra  condition 
gives  great  difficulty  in  the  numerical  simulations  but  this  has  been  surmounted  by  a  simple  method  developed  earlier  in 
our  laboratory.  The  present  analysis  has  also  incorporated  a  certain  kind  of  imperfectness  of  the  interface  in  the  boundary 
conditions  by  the  introduction  of  a  simple  parameter,  called  the  imperfectness  parameter  here,  first  proposed  by  Wortberg  and 
his  colleague.  The  results  obtained  describe  appropriately  the  development  of  the  transient  flow  fields  due  to  the  processes  of 
evaporation  and  condensation  at  the  interfaces  across  which  the  continuous  energy  flows  are  taking  place.  Some  of  the  features 
worth  to  be  mentioned  are  that  1)  a  certain  value  of  the  latent  heat  parameter  gives  the  maxima  of  the  mass  and  energy  flows. 
This  fact,  which  is  newly  found  here,  is  due  to  the  coupling  effects  of  the  latent  heat  parameter  and  the  existence  of  the 
internal  structures  of  the  condensed  phases;  2)  the  negative  temperature  gradient  phenomenon,  a  well-known  phenomenon  at 
steady  state  in  weak  evaporation  and  condensation  problems  between  the  two  condensed  phases  having  no  internal  structures, 
seems  to  be  non-existent  in  the  present  case  with  internal  structures;  3)  the  negative  mass  flow  phenomenon,  first  noticed  and 
discussed  by  Sone  and  Onishi,  seems  to  be  non-existent  at  steady  state  but  this  surely  manifests  itself  in  a  short  period  of  time 
during  the  transitional  state  of  the  flow  fields. 


INTRODUCTION 

An  extensive  studies  on  flow  problems  associated  with  evaporation  and  condensation  processes  have  been  worked  out 
so  far  because  problems  of  this  kind  are  not  only  of  theoretical  interest  but  also  of  practical  importance,  especially 
in  various  thermo-fluid  engineering  fields.  Most  of  the  analyses  done  so  far,  however,  are  confined  to  flow  systems 
in  which  the  condensed  phases  have  no  internal  structures  at  all.  The  existence  of  the  internal  structure  within  the 
condensed  phase  may  affect  drastically  the  behavior  of  the  gas  phase  and,  hence,  the  mass  and  energy  transfers. 
In  addition,  practically  important  heat  and  energy  exchange  systems  or  transfer  systems  such  as  heat  pipes  may  be 
considered  to  be  closely  associated  with  such  flow  systems  having  internal  structures. 

Recently,  Onishi  and  Yamada  [1]  have  carried  out  a  numerical  analysis  on  the  half-space  problem  based  on  the 
Boltzmann  equation  of  BGK  type  [2],  i.e.,  evaporation  and  condensation  flows  from  or  onto  the  plane  condensed 
phase  with  temperature  field  as  its  internal  structure,  revealing  the  effects  of  the  various  parameters  involved,  the 
parameters  associated  with  the  thermal  conductivity  and  thermal  diffusivity  of  the  condensed  phase  and  also  the  latent 
heat  of  vaporization,  on  the  flow  fields.  Among  others,  the  role  of  the  latent  heat  parameter  was  found  to  be  totally 
different  from  the  case  in  which  the  condensed  phase  has  no  temperature  field,  i.e.,  the  condensed  phase  has  the 
infinitely  large  thermal  conductivity.  To  be  specific,  for  a  certain  range  of  values  of  the  latent  heat  parameter,  the  phase 
change  processes  become  strong  and,  therefore,  the  associated  mass  and  energy  flows  become  large.  This  indicates 
that  there  may  exist  a  certain  value  of  the  latent  heat  parameter  which  gives  the  maxima  of  the  mass  and  energy  flows. 
This  fact  is  in  distinctive  contrast  to  what  we  know  from  the  studies  done  so  far  about  the  mass  and  energy  flows  when 
the  condensed  phases  have  no  temperature  fields  as  their  internal  structures.  In  the  latter  case,  the  mass  and  energy 
flows  become  large  with  the  values  of  the  latent  heat  parameter.  In  view  of  this,  theoretical  and  practical  interests 
will  surely  be  aroused  as  to  the  effects  of  those  parameters  mentioned  above  on  the  two-surface  problem  in  which 
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two  condensed  phases  are  involved,  i.e.,  on  the  flow  fields  confined  by  the  two  condensed  phases  having  temperature 
fields  as  their  internal  structures  (having  finite  thermal  conductivities).  This  is  what  the  present  work  is  concerned 
with.  The  reason  is  that,  in  this  two-surface  problem,  when  the  condensed  phases  have  no  internal  structures  at  all, 
the  phenomena  of  the  so-called  negative  temperature  gradient  [3] -[4]  and  the  negative  mass  How  [4]  may  appear.  The 
former  phenomenon  manifests  itself  for  values  of  the  latent  heat  parameter  larger  than  a  certain  definite  value  when 
the  phase  change  processes  are  fairly  weak.  The  latter  does  for  its  extremely  smaller  values. 

Finally,  a  certain  kind  of  imperfectness  of  the  interfaces  is  also  taken  into  account  in  the  present  study  through  the 
boundary  conditions  for  the  distribution  functions  at  the  interfaces.  This  is  expressed  in  terms  of  a  parameter,  called 
the  imperfectness  parameter  here,  first  introduced  by  Wortberg  and  his  co-workers  [5].  The  effects  of  this  parameter, 
however,  will  not  be  discussed  here  owing  to  the  limited  pages. 


FORMULATION  OF  THE  PROBLEM 


Consider  a  vapor  between  the  two  condensed  phases  with  finite  thermal  conductivity  placed  in  parallel.  The  condensed 
phases  occupy  regions  of  -D\  ^  v  ^  0  and  L  ^  v  ^  L  +  D2 ,  respectively.  Let  the  interface  surfaces  between  the  gas 
phase  and  the  condensed  phases  be  at  v  =  0  and  x  =  L.  Initially,  the  vapor  phase  and  its  condensed  phases  are  in 
complete  equilibrium  at  a  temperature  Tq.  Let  the  pressure  and  density  (or  number  density)  of  the  vapor  at  this  state 
be  Po  and  po  (or  No),  respectively.  Suppose  that  at  a  certain  time,  say  t-  0,  the  temperatures  of  the  edge  surfaces  of 
the  condensed  phases  at  x=  -D\  and  v  =  L  +  T>2  are  suddenly  changed  to  Tc\  and  Tc 2,  respectively.  The  heat  flows 
then  occur  through  the  condensed  phases  and  after  a  certain  elapse  of  time  the  temperatures  of  the  interface  surfaces 
start  to  be  changing  with  time.  This  leads  to  the  onset  of  phase  change  processes  at  the  interfaces,  giving  then  rise  to  a 
transient  motion  of  the  vapor  in  between.  This  is  what  the  present  study  is  concerned  with 

The  equation  of  heat  conduction  with  constant  substance  properties  is  assumed  to  be  able  to  govern  the  temperature 
field  within  each  of  the  condensed  phases,  which  may  be  written  as 
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where  t  is  the  time;  v  is  the  coordinate;  kc  is  the  thermal  diffusivity  of  the  condensed  phases  and  7)  (  /  =  1,2) 
represents  the  temperature  fields.  For  the  description  of  the  motions  of  the  vapor  in  the  present  problem,  on  the  other 
hand,  the  Boltzmann  equation  of  BGK  type  [2]  is  used,  which  may  be  written  as 
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where  (£*,  fy,fz)  is  the  molecular  velocity  vector;  /  is  the  molecular  velocity  distribution  function,  Fe  being  the 
local  Maxwellian  distribution  characterized  by  the  local  fluid  dynamic  quantities;  N,  u,  T ,  P  and  p  are,  respectively, 
the  number  density,  the  velocity,  the  temperature,  the  pressure  and  the  density  of  the  gas;  m  is  the  molecular  mass; 
k  is  the  Boltzmann  constant  and  R  =  k/m  the  gas  constant  per  unit  mass  of  gas.  vc  is  a  constant  associated  with  the 
collision  frequency  ( Nvc  is  the  local  collision  frequency)  and,  hence,  can  be  calculated  either  from  the  viscosity  p 
or  from  the  thermal  conductivity  A  of  the  gas  at  a  certain  reference  state,  say,  at  the  initial  equilibrium  state,  by  the 
following  relation 
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the  suffix  0  being  understood  to  indicate  the  quantities  at  the  initial  state.  It  may  be  noted  here  that,  in  the  BGK  model 
equation,  the  relation  A  =  (5/2 )Rju  holds  and,  hence,  the  Prandtl  number  Pr  is  unity  for  this  model  equation. 

The  initial  conditions  for  the  present  problem  may  be  specified  as  follows: 

f\=f2  =  T0  (-Di^jc^O,  L^x^L  +  D2)  (7) 


for  the  temperatures  of  the  condensed  phases  and 
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for  the  distribution  function  of  the  molecules  of  the  gas  phase,  the  components  of  the  molecular  velocity  vector 
(£c,  having  all  possible  values.  The  boundary  conditions  at  the  edges  of  the  two  condensed  phases  and  at 

the  interface  surfaces  between  the  gas  phase  and  the  condensed  phases  become  as  follows: 
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for  molecules  leaving  each  of  the  interface  surfaces,  respectively,  where  7Vz,  which  is  an  unknown  parameter, 
represents  the  temperature  of  the  interface  surface  itself  and  will  eventually  be  determined  as  part  of  the  solution. 
Nwi,  which  is  given  by 

Nwi  =  &cNwi  +  (l-&c)NDi  ( 0  ^  ac  ^  1 )  (12) 


is  the  number  density  for  those  molecules  emitted  from  the  interface,  ac  being  a  parameter  associated  with  a  certain 
kind  of  the  imperfectness  of  the  interface  (see  [5])  and  is  called  the  imperfectness  parameter  here.  Nwi  is  the  saturated 
vapor  number  density  at  the  temperature  Twi  of  the  surface  of  the  condensed  phase  to  be  determined  by  the  Clapeyron- 
Clausius  relation  as 


Nwi  = 


Pwi 

kTwi’ 


Pwi 


=  exp 


(13) 


where  T  is  a  non-dimensional  parameter  defined  by  T  =  HlURTo),  Hl  being  the  latent  heat  of  vaporization  per  unit 
mass.  The  other  quantity  Not ,  on  the  other  hand,  is  the  number  density  for  molecules  reflected  from  the  condensed 
phase  determined  by  the  condition  of  no  net  mass  Row  across  the  phase  boundary.  Nuu  therefore,  is  obtained  from  the 
distribution  function  for  incident  molecules  to  the  condensed  phase  as 
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the  upper  sign  applying  to  the  interface  at  x  =  0  and  the  lower  to  that  at  x  =  L,  respectively.  It  may  be  noted  that 
Not  is  also  an  unknown  constant  to  be  determined  as  part  of  the  solution.  In  addition  to  the  above  specification  of  the 
conditions  at  the  boundary  surfaces,  the  condition  of  the  continuity  of  energy  flow  across  each  of  the  interfaces  has  to 
be  imposed,  which  may  be  written  as 

+  Ac^-  =  ±E±pu(hL-h)  (15) 

ox 

where  Ac  is  the  thermal  conductivity  of  the  condensed  phase  itself  (assumed  constant),  h  is  the  enthalpy  per  unit  mass 
of  the  gas  and  h-  cpT,  cp  being  the  specific  heat  at  constant  pressure.  E  is  the  component  of  the  energy  flux  of  the 
gas  in  the  x -direction  defined  by 
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which  is  here  to  be  evaluated  at  the  interfaces. 
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CHARACTERISTIC  PARAMETERS 


For  the  analysis  of  the  present  problem,  we  have  introduced  L  as  the  length  scale  and  to  as  the  time  scale  taken  as 

-  =(Z)1/2  L=  1  =  n7x 

T° "  (2RTo)1/2  \2)  co  N0vc  P0  ^  ' 

where  co,  taken  as  the  velocity  scale  here,  is  the  sound  speed  of  the  gas  at  the  initial  state  defined  by  co  =  (yRTo)i/2, 
y  being  the  specific  heat  ratio  (y  =  5/3  here).  With  the  fluid  dynamic  quantities  at  the  initial  state  together  with 
these  length,  velocity  and  time  scales,  the  system  of  the  governing  heat  conduction  and  kinetic  equations  and  the 
initial  and  boundary  conditions  is  appropriately  nondimensionalized,  giving  the  following  non-dimensional  parameters 
characterizing  the  present  flow  field 
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where  do  and  kq  =  Ajq/(pocp)  are  the  thermal  conductivity  and  the  thermal  diffusivity  of  the  gas  at  the  initial  state, 
respectively.  Kn  is  the  Knudsen  number  defined  by  Kn  =  Iq/L,  Iq  being  the  mean  free  path  of  the  gas  molecules 
at  the  initial  state  defined  by  Iq  =  (po/Po)(8RTo/n)1/2.  It  may  be  noted  that  the  Reynolds  number  Re  defined  by 
Re  =  poCqL/jjo  is  related  here  to  the  Knudsen  number  as  Re  =  (8y/7r)1/2(l/ Kn). 


RESULTS 

The  system  of  the  governing  kinetic  equations  has  been  solved  numerically  subject  to  the  initial  and  boundary 
conditions  with  a  simple  difference  scheme  applied.  Here  only  a  part  of  the  results  will  be  shown.  A  sample  of  the 
development  of  a  transient  flow  field  to  its  steady  state  is  shown  in  Fig.  1  and  2,  where  the  temperature  jumps  at  the 
interfaces  x  =  0  and  x  =  L,  respectively,  can  clearly  be  visible  at  the  final  state  (t/ro  =  120.0  is  considered  to  be  the 
final  steady  state). 


FIGURE  1.  Transient  distributions  of  the  pressure  and  the  density.  The  interface  surfaces  are  at  x/L  =  0  and  x/L  -  1,  respectively. 
D\/L  =  D2/L  =  0.1,  Tci/Tq  =  2.0,  Tc2/Tq  =  1.0,  T  =  3,  Ac/A0  =  80.0,  kc/k0  =  0.50,  Kn  =  0.005  (Re  =  412.03),  crc  =  1.0.  The  numbers 
in  the  figures  indicate  the  time  t/ro. 

The  effects  of  the  latent  heat  parameter  T  coupled  with  the  thermal  conductivity  ratio  Ac/Aq  on  the  flows  are 
virtually  the  same  as  in  the  half-space  problem  done  by  Onishi  &  Yamada  [1],  i.e.,  T  — >  small  or  Ac/Ao  — >  large, 
the  temperature  gradient  &T /dx  within  the  condensed  phase  can  become  small,  which  means  the  possible  existence 
of  fairly  large  temperature  jump  at  the  interface  leading  to  the  fairly  strong  phase  change  process  taking  place.  The 
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FIGURE  2.  Transient  distributions  of  the  temperature  and  the  velocity  (continued  on  Fig.  1).  The  temperature  profiles  within  the 
condensed  phases  are  at  t/ro  =  120.0 


mass  and  energy  flows  naturally  become  large  as  T  — >  small.  However,  this  applies  down  to  a  certain  definite  value 
of  T ;  when  T  becomes  smaller  than  this  definite  value,  the  mass  and  energy  flows  begin  to  decrease  again.  Figure  3 
shows  this  behavior,  indicating  that  the  existence  of  the  maxima  of  the  mass  and  energy  flows  at  a  certain  value  of  the 
parameter  T.  This  will  be  a  new  feature  among  the  roles  of  the  latent  heat  parameter. 

Another  feature  of  the  present  problem  is  that  the  negative  mass  flow  manifests  itself  during  a  certain  period  of  time 
when  the  latent  heat  parameter  T  becomes  extremely  small  but  it  disappears  soon  as  time  goes  on.  A  typical  example 
of  this  is  shown  in  the  graph  on  the  right  hand  side  of  Fig.  4,  where  the  transient  mass  flow  distributions  are  also  listed. 
The  existence  of  this  negative  mass  flow  phenomenon,  which  has  been  noticed  at  the  steady  state  and  discussed  by 
Sone  &  Onishi  [4],  may  be  explained  as  follows:  consider  the  difference  between  the  mass  fluxes  due  to  the  molecules 
leaving  the  interfaces  at  x  =  0  and  x  =  L,  respectively,  which  is  in  this  case 

_  1  Po  (  Pwi  Pw2 

pu\x=0  pu\x=L  -  ^  {2RTo)1/2  [ 

because  the  saturated  vapor  pressures  Pwi  and  Pwi  remain  nearly  equal  even  when  tw2±Pw\  for  very  small  values 
of  the  latent  heat  parameter  T  (see  Eq.  (13)).  Molecular  collisions,  of  course,  affect  this  difference  in  the  mass  flow. 
However,  the  above  expression  can  be  applied  locally  near  the  interface,  in  which  case  the  interface  at  x-L  above,  for 
example,  becomes  the  position  near  the  interface  at  x  -  0 ,  say,  somewhere  outside  the  Knudsen  layer  (temperature  of 
the  vapor  near  the  interface  at  x  =  0,  of  course,  is  lower  than  7Vi  in  the  case  considered). 

Finally,  as  to  the  negative  temperature  gradient  [3] -[4],  a  well-known  phenomenon  for  two-surface  problems,  seems 
not  to  manifest  itself  in  the  present  cases  with  internal  structures. 


negative  for  Twi  >  Twi  for  Pwi  —  Pwi  (19) 
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FIGURE  3.  The  steady  state  mass  flow  pu/poco  and  energy  flow  E*  =  Ej P^dRT^'1  versus  the  latent  heat  parameter  T. 
D\/L  =  D2/L  =  0.1,  rci/7o  =  2.0,  Tc2/T0  =  1.0,  Kn  =  0.005  (Re  =  412.03),  ac  =  1.0.  The  white  circle:  Ac/A0  =  80.0,  kc/k0  =  0.50; 
The  balck  circle:  AJ  A o  =  30.0,  kc/kq  =  0.32.  The  existence  of  the  maxima  of  the  mass  and  energy  flows  can  be  expected. 


FIGURE  4.  A  transient  to  steady  state  mass  flow  distributions.  D\/L  -  D2/L  =  0.1,  Tci/Tq  =  2.0,  TC2ITq  -  1.0,  T  =  0.1, 
Ac/Ao  =  80.0,  kc/kq  =  0.5,  Kn  =  0.005  (Re  =  412.03),  ac  =  1.0.  The  numbers  in  the  figures  indicate  the  time  t/r 0.  The  graph 
on  the  right,  which  is  the  enlarged  portion  of  the  region  near  the  interface  at  x/L  =  0,  indicates  the  negative  mass  flow  persisting  up 
to  about  t/ro  =  4.5. 
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